Swift and XMM-Newton Observations of the Extraordinary 
GRB 060729: More than 125 days of X-ray afterglow 



Dirk Grupe^, 
grupeOastro . psu . edu 

Caryl Gronwall^, Xiang-Yu Wang^'^, Peter W.A. Roming^, Jay Cummings^, Bing Zhang^, 
Peter Meszaros Maria Diaz Trigo^, Paul T. O'Brien^, Kim L. Page^, Andy Beardmore^, 

Olivier Godet^, Daniel E. vanden Berk^, Peter J. Brown^, Scott Koch^, David Morris^, 
Michael Stroh^, David N. Burrows^, John A. Nousek^, Margaret McMath Chester^, Stefan 
Immler*'^, Vanessa Mangano^°, Patrizia Romano^^, Guido Ghincarini^^ , Julian Osborne^, 

Takanori Sakamoto^, Neil Gehrels^ 



ABSTRACT 

We report the results of the Swift and XMM-Newton observations of the S'wi/i-discovered 
long Gamma-Ray Burst GRB 060729 (Tgo=115s). The afterglow of this burst was exceptionally 
bright in X-rays as well as at UV/Optical wavelengths showing an unusually long slow decay 
phase (a=0.14±0.02) suggesting a larger energy injection phase at early times than in other 
bursts. The X-ray light curve displays a break at about 60 ks after the burst. The X-ray decay 
slope after the break is Q!=1.29±0.03. Up to 125 days after the burst we do not detect a jet break, 
suggesting that the jet opening angle is larger than 28°. In the first 2 minutes after the burst (rest 
frame) the X-ray spectrum of the burst changed dramatically from a hard X-ray spectrum to a 
very soft one. We find that the X-ray spectra at this early phase can all be fitted by an absorbed 
single power law model or alternatively by a blackbody plus power law model. The power law fits 
show that the X-ray spectrum becomes steeper while the absorption column density decreases. In 
the blackbody model we find that the temperature changes from kT = 0.6 keV at 85 s after the 
burst to 0.1 keV at 160 s after the burst in the rest frame. In Swift^s UV/Optical telescope the 
afterglow was clearly detected up to 9 days after the burst in all 6 filters and even longer in some 
of the UV filters with the latest detection in the UVWl 31 days after the burst, which is one of 
the latest detections of an afterglow by the UVOT. A break at about 50 ks is clearly detected in 
all 6 UVOT filters from a shallow decay slope of about 0.3 and a steeper decay slope of 1.3. The 
spectral energy distributions show that there is no spectral development between X-rays and the 
UV/optical throughout the observation. In addition to the Swift observations we also present 
and discuss the results from a 61 ks Target-of-Opportunity observation by XMM-Newton starting 
about 12 hours after the burst. These observations show a typical afterglow X-ray spectrum with 
/3x=l-l and absorption column density of 1 x 10^^ cm^^. 

Subject headings: GRBs: individual(GRB 060729) 
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1. Introduction 

Gamma-Ray Bursts (GRBs) are the most pow- 
erful explosions in the present-day Universe. With 
the launch of the Swift Gamma- R ay Burst Ex- 
plorer Mission ( Gehrels et al.ll2003 ) in November 
2004 a new era in GRB science has started. Swift 
is able to observe the afterglow of a burst with 
its nar row-field instruments, the X-ray Telescope 
(XRT, iBurrows et allbOOSaf) and the UV/Qptical 
telescope (UVOT: iRoming et al.ll2005l ). typically 
within 2 minutes after th e detection by the Burs t 



Alert Telescope (BAT, iBarthelmv et al.l l2005[ ) 



Several phenomena were discovered by Swift, like 
the occurrence of giant fl ares during the first 
1000 s after the bu rst (e.g. iBurrows et al.|[2005bt 
Falcone et al.ll2006 ) or the canonical Hght curves of 



GRB afterglows (jNousek et al.l 120061 : IZhang et al 
20061 ). 



GRB 060729 was discove red by Swift on 2006 
July 29 ( Grupe et al. 20061 ) as one of the bright- 
est bursts ever detected in X-rays by Swift 
(With the brightest GRB in XRT so far be- 
ing GRB 0611 21: iPage et all |2007|) . Besides 



GRB050525A (iBlustin et al.ll2005D . GRB 060218 
(ICampana et al.l2006l) . g"rB 060614 (iMangano et al 



in preparation! ) ■ and GRB 061121 ( Page et al 
20071), GRB060729 is the burst with the best 



UVOT follow up even up to 9 days after the burst 
in all 6 UVOT filters and even longer in some of 
the UV filters. In the XRT the afterglow has been 
detected by more than 125 days after the burst. 
This is the longest follow-up observation with a 
detection of an afterglow ever performed by Swift. 
Similar coverage, but shorter, has only been per- 
formed for GRBs 050416A, 060319, and 060614. 
Even though the burst was bright in the BAT, 
it was too faint to be detected by Konus-Wind 
(Frederiks, priv. comm.). 
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Even though the sun angle was small in RA 
(2.2 h), due to the declination of Dec=-62° 
the afterglow was circumpolar for most south- 
ern observatories, and was observed by the ESO 
VLT using F0RS2 and by G emini South us- 



ing GMOS (|Thoene et al. 

z= 0.54 was determined from the optical spectra 



20061 ). A redshift of 



bv iThoene et all l|2006h . GRB 060729 was also 
observed by ROTSE Ilia, located at the Siding 
Spring Observatory, Australia, by lOuimbv et al. 
( 20061 ). who reported an initial upper limit of 16.6 



mag 64 s after the BAT trigger. They were able 
to detect the afterglow with ROTSE Ilia up to 
175 ks after the burst when it was still at 19.4 



mag (jQuimbv fc Rvkoff 1 120061) . decaying with a 
slope a=0.23. Cobb fc Bailvn ( 20061 ) measured 
a decay slope in the I band of ai=1.5 based on 
GTIO 1.3m SMARTS observations between 4.6 to 
17.6 days after the burst. 

The paper is organized as follows: In |2]we de- 
scribe the observations and the data reduction. 
In fj3] we present the data analysis. The discus- 
sion of our results is given in 21 Throughout 
the paper decay and energy spectral indices a 
and /3 are defined by F,y{t,iy) oc {t ~ to)~°'iy~^, 
with to the trigger time of the burst. Luminosi- 
ties are calculated assuming a ACDM cosmology 
with f7M=0.27, flA=0.73 and a Hubble constant 
of Hq—71 km s~^ M pc~^ us i ng th e luminosity 
distances Dl given bv iHogg I (1 19991 ) resulting in 
£'l=3120 Mpc. All errors are la unless stated 
otherwise. 

2. Observations and data reduction 

The Swift BAT triggered on the pre-cursor 
of GRB 060729 at 19:12:29 UT o n 2006 July 29 
()Grupe et al.l l2006t iGrupe I l2006f ) . Swift's X-ray 
Telescope began observing the afterglow 124 s af- 
ter the trigger. The UVOT started the observa- 
tions 135 s after the BAT trigger. 

The Swift XRT observed GRB 060729 in the 
Windowed Timing ( WT) and Photo n Gounting 
(PC) observing modes (jHiU et al.ll2004l ). The XRT 
data were reduced by the xrtpipeline task version 
0.10.4. The WT mode data at the beginning of 
the XRT observation had to be treated with spe- 
cial care. During the first 20s after the start (130- 
150s after the BAT trigger) of the WT XRT ob- 
servation the satellite was still settling causing the 
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target to move on the XRT C CD towards the dea d 
columns at DETX=319-321 (|Abbev et al.ll2006h . 
In order to correct for the photon losses due to 
these dead columns we measured the offset of the 
source at each second and calculated a correction 
factor according to the losses of the WT mode 
Point Spread Function (PSF). After 150s after the 
trigger all WT mode data were corrected by the 
same factor. Source and Background photons were 
selected by XSELECT version 2.4 in boxes with a 
length of 40 pixel. For count rates > 150 counts 
s~^, however, in the WT mode the data had to 
be corrected for pileup. In order to correct for 
pileup in the light curve and for the spectral anal- 
ysis and the determination of the hardness ratic{3 
we excluded the central regions of the PSF at the 
source po sition depending on th e count rate as de- 
scribed in I Romano et al. ( 2006f l. For the PC mode 
data the source photons were selected in a circu- 
lar region with a radius of r=59 and the back- 
ground photons in a circular region close by with 
a radius r=176 in the first segments. For the 
later data the radii were reduced to 47 and 24 
for the source, and 137 and 96 for the back- 
ground. For the spectral data only events with 
grades 0-2 and 0-12 were selected with XSELECT 
for the WT and PC mode data, respectively. Note 
that the source photons for the spectral analysis 
of the PC mode data of the first orbit were se- 
lected in a ring with an inner radius of 16'.'5 and 
an outer radius of 71"0 in order to avoid the effects 
of pil e-up (e.g. Pagani et al. 20061 : Vaughan et al 
20061 ). The spectral data were re-binned by us- 



ing grppha version 3.0.0 having 20 photons per 
bin. The spe ctra were ana. lyzed with XSPEC 
version 12.3.0 (|Arnaudlll996h . The auxiliary re- 
sponse files were created by xrtmkarf and corrected 
using the exposure maps, and the standard re- 
sponse matrices swxwt0to2-20010101v008.rmf and 
swxpcOto 12-2001 01 OlvOOS.rmf were used for the 
WT and PC mode data, respectively. All spectral 
fits were performed in the observed 0.3-10.0 keV 
energy band. For the errors of the spectral fit 
parameters we used the standard 
in XSPEC, that is equivalent to a 90% con- 
fidence region for a single parameter. 



^The hardness ratio is defined by HR = (hard — 
soft) /{hard + soft) where soft is the counts in the 0.3- 
1.0 keV band and hard is the counts in the 1.0-10.0 keV 
band, respectively. 



Background-subtracted X-ray flux light curves 
in the 0.3-10.0 keV energy range of the Swift ob- 
servations were constructed using the ESQ Mu- 
nich Image Data Analysis Software MIDAS (ver- 
sion 04Sep) and by an IDL program that corrects 
for PSF losses, in particular when the source is 
located on one of the dead columns on the XRT 
CCD detector. The light curve was binned as fol- 
lows: the WT mode data with 1000 photons per 
bin, and the pc mode data with 200 counts per 
bin in the first days after the burst and 20 or 10 
at the end of the observations. The count rates 
were converted into unabsorbed flux units using 
energy conversion factors (ECF) which were de- 
termined by calculating the count rates and the 
unabsorbed fluxes in the 0.3- 10.0 keV energy ban d 
using XSPEC as described in lNousek et al.r(|2006f ). 
The XRT data at the beginning of the observa- 
tion show dramatic spectral changes and require 
specific ECFs for each time bin. The later data, 
however, do agree with a typical afterglow spec- 
trum and the count rates were converted by one 
ECF=5xlO~^-^ ergs s~-^ cm~^ (counts s~^)~-^. 

The Swift UVOT observations of GRB060729 
began with the automated "GRB" sequence, 
which provided finding chart images in White (100 
s) and V (400 s), then began cycling through all 
6 UV and optical filters starting 739 s after the 
trigger. The source data of these early White and 
V images were extracted from a circle with a ra- 
dius of 6 . GRB 060729 remained detectable in 
all 6 filters for more than 9 days after the trig- 
ger, then for 12 days after the trigger in UVW2 
(Ac — 1930A) and for 31 days after the trigger 
in UVWl (Ac = 2510A). During the first days 
after the burst each observation in each single or- 
bit was analyzed. For the later data the images 
were coadded with uvotimsum in order to improve 
the signal to noise ratio. The data were analyzed 
with the UVOT software tool uvotsource. Due 
to the bright F3V star HD 45187 (9.4 mag in B, 
107'.'5 away from GRB 060729) extra caution had 
to be taken for the source extraction and the back- 
ground subtraction. We chose a selection radius 
of 4 for the source and 8 for the background in 
all filters placed at a position nearby on the rim of 
the bright star's halo as displayed in Figurc[TJ In 
order to correct for losses due to this small source 
extraction radius, we did an aperture correction 
with V=0.03 mag, B=0.05 mag, U=0.05 mag. 
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UVW1=0.17 mag, UVM2(Ac = 2170A)=0.15 
mag, and UVW2=0.15 mag. All values plotted 
and listed in this paper take these corrections into 
account. The data, however, are not corrected for 
Galactic reddening which is i5B-v=0.050 mag 
(ISchleeel et al.lll998l) in the direction of the burst. 

GRB 060729 was also observed by XM M-Newton 
(IJansen et al.l 200ll) for a total of 61 ks ( Schartel 
20061 : ICampana k De Luca Il2006f ) . XMM-Newton 
started observing the afterglow of GRB 060729 on 
2006 July 30 07:41 (44.9 ks after the trigger) to 
2006 July 31 01:04 UT (107.5 ks after the trigger). 
In t he European Photon Imaging Camera (EPIC) 
pn (jStriider et al.l 120011 ) the total observing time 
was 59.6 ks using the medium light blocking filter. 
However, due to high particle background during 
part of the observation only 42.3 ks were used 



The observations in the EPIC MOS (Turner et al 



20011 ) were for a total observing time of 61.2 ks. 
The MOSl was using the medium filter while the 
M0S2 observations were performed with the thin 
filter. The total observing time in the Reflection 
Grating Spectrometers (RGS; Iden Herder et al.l 
2001.) was 615 ks In the Optical Monitor (OM; 
Mason et al.ll2001l) the afterglow was observed for 
8.3 ks in U, 5x 4 ks in UVWl, 3x 4 ks in UVM2, 
and 2x 4 ks with the optical grism. Note, that 
the OM UVWl and UVM2 filters are slightly dif- 
ferent compared with the UVOT UV filters. The 
XMM-Newton data were reduced with the latest 
SAS version xmmsas.20060628.1801-7.0.0. 

3. Data Analysis 

3.1. Position of the Afterglow 

The position of the afterglow measured from 
the UVOT UVWl coadded image is 
RA (J2000) = 06^21™31';86, 
Dec (J2000) = -62°22'12'.'5 

with a 1 error. This position is consistent with 
the initi al analysis of the White and V filter 



analysis ( Immler 2006( ). The UVWl position is 



I'.'l away from the X-ray position RA (J2000) = 
06^21™31!75, Dec (J2000) = -62°22'13'.'3 (with 
a 3'.' 5 90% confidence error), which was measured 
for the XRT PC mode data of segment 001 using 
the new t eldef file swx20 m02v 001.teldef as de- 
scribed in iBurrows et al.l (|2006al ). This position 
deviates by 3" 2 from the refined position given in 
Grupe ( 2006f l. The most likely reason for this 



differen c e is t hat for the X-ray position given in 
Grupe ( 20061 ) only the PC mode data of the first 
orbit were used. This is due to the burst dur- 
ing the first orbit being placed on one of the bad 
columns on the XRT CCD which makes the de- 
termination of a position difficult. Figure[T] dis- 
plays the UVWl image of the field of GRB 060729. 
The circle in the upper right inserted image is the 
3'.' 5 XRT error radius of the X-ray position given 
above. 

3.2. BAT data 

Figured] displays the BAT light curves in 
the 15-25 keV, 25-50 keV, 50-100 keV, and 
15- 100 keV bands (top to bottom) with To= 
2006 July 29 19:12:29 UT (Spacecraft clock 

175893150.592 )^ GRB 060729 had a r9o=115s 

2006f l. Partly the Tgo is so long be- 



( Parsons et al. 



cause the trigger was on the pre-cursor. After the 
initial first peak (pre-cursor) which was detected 
by the BAT and which triggered the observation 
the burst drops back down to the background 
level. However, two giant peaks are observed at 
about 60s after the trigger, of which the first one 
is harder than the second. There is a third peak 
which peaks about 120s after the trigger. This 
is the peak of which we see the end of the decay 
in the XRT (Figure [3|) . For the spectral analysis 
the BAT data were divided into 5 bins as listed 
in Tabled] The first peak is the initial peak the 
BAT triggered on GRB 060729. As shown in Ta- 
blejl] the following two peaks that occur between 
70 and 124s after the burst are by a factor of 3 
stronger than the initial peak. The two last peaks 
(124- 190s after the burst, XRT fiare 1 and 2 in 
Table]!]) are also observed simultaneously in the 
XRT. These data will be discussed in the XRT 
section. 

TableU] lists the results of the spectral analy- 
sis of the 5 peaks. All spectra were fitted by a 
single power law model. The initial peak has a 
hard spectrum with a 15-150 keV energy spec- 
tra index /3i5_i5okoV=l-05lo.32- The two strong 
peaks between 70-124s after the burst show an 
interesting spectral behavior. While the first of 
these peaks (70-88s after the burst) has a rather 
hard spectral slope with /3i5_i50kcV=0.59±0.11, 
the second of these peaks (88- 124s) was softer with 
/3i5-i50keV=0.90±0.11. The total fluence in the 
observed 15-150 keV band is 2.7x10^^ ergs cm"^ 
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( Parsons et al ] |2006h and 7.2x10 ergs cm and 
1.7x10"^ ergs cm"^ in the rest-frame 1 keV - 1 
MeV and 1 keV - 10 MeV bands, respectively, 
adding all BAT spectra together (Tabled]) and as- 
suming the same power law spectrum as in the 
15-150 keV band without any break. With the 
redshift z=0.54 this converts into an isotropic en- 
ergy in the rest-frame 1 keV - 1 MeV and 1 keV 
- 10 MeV band of Ei^o = 6.7 x lO^^^ ergs and 
i^iso = 1.6 X 10^^ ergs, respectively. Because we 
lack observations of the break energy inbreak and 
the 7-ray spectrum at higher energies by Konus- 
Wind the 1 keV - 10 MeV band Eiso value is an 
upper limit of the true isotropic energy. 

3.3. X-ray data 

3.3.1. Temporal analysis 

FigureO shows the combined BAT and XRT 
light curve. The light curve clearly shows that 
XRT began observing the GRB at the beginning 
of the fourth peak seen in the BAT. The com- 
bined BAT -|- XRT light curve wa s constructed 



as described in O'Brien et al. I (|2006h . Due to the 



dramatic spectral change within the three min- 
utes of the WT observation, we applied an ECF 
for each individual bin assuming power law model 
corrected for absorption with the parameters listed 
in Table[21 However, for the BAT data we applied 
only one ECF that reflects the main spectrum. 

Figure|4] displays the Swift XRT light curve, 
with WT mode data as triangles and PC mode 
data as crosses. The vertical dashed lines in 
the figure mark the start and end times of the 
XMM-Newton observations. The general behavior 
of the light curve can be described as follows: af- 
ter the initial steep decay with a decay slope ai 
= 5.11±0.22 the light curve flattens at Tbrcak,! = 
530s±25s with a decay slope a2 = 0.14±0.02. At 
Ibreak,2 = 56.8 ksiblOks the light curve of the af- 
terglow breaks again and continues decaying with 
a decay slope as = 1.29±0.03. The definitions of 
the decay slopes fo llow the descriptions giv en in 
Nousek et all (|2006l ) and lZhang etaD (|2006f ). We 
do not detect a jet break even 125 days after the 
burst. The last Str detection of the X-ray after- 
glow was obtained between 2006 November 21 to 
December 01 with a total exposure time of 69.9 
ks. The afterglow was still observed by Swift un- 
til December 27 for a total of 63.5 ks. However 



these observations were interrupted by several new 
bursts and at the end only a 3cr upper limit of 
2.1x10"^^ ergs s~^ cm~^ could be obtained. It 
was dropped from the Swift schedule after 2006 
December 27 because it was not detectable any- 
more with the XRT within a reasonable amount 
of observing time. 

3.3.2. Spectral Analysis 

Dramatic spectral change at the beginning 

In order to examine the spectral behavior in more 
detail, source and background spectra were cre- 
ated for each bin, except for the first 10 bins 
(T - ro=130~150 s). The WT mode data were 
divided into 21 bins with 1000 source photons in 
each bin. Because of the high count rate at the be- 
ginning of the observations we applie d the method 
as described in Romano et al. ( 20061 ) to avoid the 
effects of pileup. However, this procedure reduced 
significantly the number of source photons in each 
single spectrum. We, therefore, combined two bins 
into one for bins 1-^2, 3-h4, 5-h6, 7+8, and 9+10 
to increase the signal to noise ratio. 

Each of the 16 spectra were fitted by an ab- 
sorbed single power law, blackbody plus power 
law, and a power law with exponential cutoff 
model. The results of these spectral fits are listed 
in Table [21 Figure[5] shows plots of the WT mode 
data; from top to bottom the count rate, hard- 
ness ratio, the X-ray spectral slope /3x from an 
absorbed single powerlaw fit with a free fit absorp- 
tion column density iV/f, blackbody temperature 
kT (in keV) from the blackbody plus power law 
fit, the blackbody radiu^ i?bb and the break en- 
ergy £^brcak of a powcr law with exponential cutoff 
model. The hardness ratio changes from HR=0.6 
at the beginning of the observation to HR=— 0.56 
at the end, indicating a dramatic evolution in the 
X-ray spectrum within 2 minutes of observing time 
which translates into 1.3 minutes in the rest frame. 

While the spectra during the first 20 s of the 
WT mode observation are well-fitted by a single 
power law model, after about 150s after the burst 
the data are better fit by a blackbody plus power 
law spectrum. For the absorbed power law fits. 



^The blackbody radii were derived for each bin from the 
relation L ■ 
constant 



: A-kR^^ctT'^, where a is the Stcfan-Boltzmann- 
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all parameters were left free. We found that while 
the spectra at the beginning of the observation 
were rather hard with /?x=l-5 and A'h=4 x 10^^ 
cm"^, the spectra became very soft with (3x^3.0 
and A'^H=1 x 10^^ cm~^. Note that the absorp- 
tion column density TVh decreases while the energy 
spectral index /3x becomes steeper - the opposite 
than what is expected if the spectral slope and 
the absorption column density were just linked in 
the fitting program. However, note that especially 
during the later bins, the spectra are not well fit 
by a single power law and do require more com- 
plicated models. 

For the blackbody plus power law model, the 
absorption parameter was fixed to the Galactic 
value (4.82x10^0 cm"^ iDickev fc Lockma^ll990l ) 
and the hard energy spectral slope to /3x=l-0. The 
blackbody temperature changes dramatically from 
fcT=0.56 kcV at the beginning of the XRT WT ob- 
servation to 0.11 keV at the end, accompanied by 
an increase of the blackbody radius from 2.5 x 10^^ 
cm to 16x10^^ cm. Fitting the data with an ab- 
sorbed blackbody plus power law model with the 
absorption column density at z=0 set to the Galac- 
tic value and at z~Q.54 set to 1 x 10^^ cm~^ (see 
the discussion about the XMM-Newton spectral 
analysis) results in similar values for the temper- 
ature. The only differences are that the tempera- 
tures tend to be lower by 40 cV and the normal- 
izations arc higher. 

The prompt emiss ion of GRBs is o ften fitted 
by a Band function ijBand et al. Il993h . We also 



tried a power law model with exponential cutoff, 
which is a surrogate for the Band model, which has 
the advantage of using fewer parameters than the 
Band model. In order to obtain better constraints 
we fixed the absorption column to the Galactic 
value. As listed in TableO typically the power law 
model with exponential cutoff does not show im- 
provements over the single power law or the black- 
body plus power law models. 

The change in the X-ray spectra is also dis- 
played in Figurc[6] which shows the spectra of bins 
# 1, 12, 14, and 21. Bin 12 is the bin before the 
small flare at 170s after the burst and bin 14 is the 
peak of that flare. 



Later PC mode data 

All PC mode data can be fitted by a power law 
model with a energy spectral slope /3x=l-2 and 
an absorption column density of about 1.5x10^^ 
cm^^. This absorption column density is signif- 
icantly above the Galactic value . The intrinsic 
absorption column density at the redshift z=0.54 
is 1.9±0.4xl02i cm-2. TablcE] lists the XRT PC 
mode observations at 20-40 ks after the burst and 
at 200 ks after the burst, so before and after the 
XMM observation. The fits to these data suggest 
no significant spectral variability before and after 
the break in the X-ray light curve around 60 ks 
after the burst. 

The XMM-Newton observations 

The combined spectra of the XMM-Newton EPIC 
pn and MOS and Swift XRT data are shown in 
Figure[7l The XRT data were selected between 
44900 s and 107500 s after the burst. The details 
of the spectral fits to these data are summarized 
in TableO At these late times, the X-ray spec- 
tra were well fitted by absorbed single power law 
models. However, as a check the spectra were fit- 
ted also by a blackbody plus power law model, but 
at these late times the power law component dom- 
inates the spectra. Therefore, we will only discuss 
the absorbed power law model fits as listed in Ta- 
bleO 

The obvious difference between the Swift XRT 
and XMM-Newton pn and MOS data is the much 
higher value of the absorption column density. 
From the free fit absorption column density at 
z=0 we measured an absorption column density 
iVH=15.7xl02° cm-2 in the SwiftXRT data. This 
value is about twice as high as what is measured 
from the XMM-Newton EPIC pn and MOS spec- 
tra. The EPIC pn is well-calibrated t o energies 
below 0.2 keV (e.g. iHaberl et aLll2003l ). We also 
applied the gain fit model within XSPEC, but it 
did not remove the discrepancy. This discrepancy 
maybe due to problems with the Swift XRT bias 
maps during the time period 2006 July 21 and 
2006 August 3. This bias map problem causes an 
offset in the gain and therefore compromises the 
spectral analysis of Swift XRT PC mode data dur- 
ing that time period. However, this gain shift does 
not affect the early Swift XRT WT mode data. 
Due to the better response of the EPIC pn at lower 
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energies we consider the absorption column densi- 
ties measured by the EPIC pn the most reliable. 
With the redshift of the burst at z=0.54 we can 
also use the X-ray spectra to determine the intrin- 
sic absorption column density at the location of 
the afterglow. The intrinsic column densities of 
all fits are in the order of 1 x 10^^ cm~^, except 
for the Swift XRT data which again show an ab- 
sorption column density about twice as high. As 
we will show later in section l531 the absorption col- 
umn density of 1 x 10^^ cm~^ is in good agreement 
to what can be derived from the spectral energy 
distribution of the afterglow. 

In addition to the EPIC pn and MOS data we 
also analyzed the 2 RGS spectra. Wc found that 
the analysis of the RGS continuum spectra agrees 
within the errors with the pn and MOS data. We 
did not find any obvious emission or absorption 
features in the RGS spectra. 

3.4. UV/Optical data analysis 

The magnitudes resulting from the UVOT data 
analysis are listed in TablejH Figure[8] shows 
the results of the UVOT photometry in compar- 
ison with the XRT. UVOT was able to follow 
the afterglow in all six filters up to 9 days af- 
ter the burst. In UVWl the afterglow was fol- 
lowed up 31 days after the burst which trans- 
lates into 20 days in the rest-frame. This is one 
of the longest intervals Swiffs UVOT has ever 
detected an a fterglow in the optic al/UV. Only 
GRBs 060218 (ICampana et al.l l2006h and 060614 
( Mangano et all in preparation! ) were detected at 
slightly later observed times than GRB 060729. 

In all bands, XRT as well as in all 6 UVOT fil- 
ters, a significant break occurs in the light curve. 
TableO lists the decay slopes a2 and as before and 
after the break time Tbrcak- Within the errors all 
break times seem to occur at about 50 ks after the 
burst (33 ks in the rest-frame), with the earliest 
break in B at about 30ks and the later breaks at 
shorter wavelengths. However, considering the un- 
certainties in the decay slopes, this is all consistent 
with an achromatic break. Note that due to a re- 
brightening of the afterglow at about 20 ks seen in 
X-rays and all 6 UVOT filters the determination 
of a2 is rather uncertain. In B the afterglow de- 
cays the slowest with a3=0.98. The decay slopes 
at shorter wavelengths are steeper with as ?a=1.3. 
Note that the flatter slope in the B-Filter is caused 



by a re-brightening at about 200 ks that is not 
seen in the other filters. By limiting the analy- 
sis to data only up to 200ks the decay slope is 
a3=1.17±0.16 which is consistent with the decay 
slopes seen in the other filters. The rebrighten- 
ing in B at about 200 ks after the burst seems 
to be real. We checked for any strong variabil- 
ity in the background but could not detected any 
strong background variation at that time. The de- 
cay slopes a2 and a.^ are consistent with the d ecay 
slopes reported by Ouimbv fc Rvkoff ( 2006f ) and 



Cobb fc Bailv^ (l2006l) 



Figure[9] displays the UVOT White and V and 
XRT light curves of the first orbit. The UVOT 
data of this period are listed in TablelHl The left 
panel of Figurc[9] displays the UVOT White filter 
event mode and XRT WT mode data. The UVOT 
White filter data were grouped into 10 s bins. The 
first UVOT White points show a decay similar to 
the XRT WT light curve. However after these few 
points the UVOT White light curve fiattens, which 
agrees with the fiare seen in the XRT data at 170s 
after the burst. At about 200s after the burst the 
afterglow starts to become brighter in the UVOT 
White, while it is still decaying in X-rays. 

The right panel of Figure[9] shows the UVOT V 
event mode data and the XRT PC mode data of 
the first orbit. The UVOT V data were grouped 
in 25 s bins and the XRT PC mode data with 25 
source photons per bin. The UVOT V light curve 
shows the afterglow fairly constant at about 17.5 
mag while the XRT PC mode light curves shows 
an initial decay until about 600s and fiattens after 
that. 

The XMM-Newton Optical Monitor observa- 
tions are summarized in Table[71 While the U filter 
results agree well with the UVOT U data as listed 
in TableSl there is a discrepancy in the UVWl 
and UVM2 filters. There maybe two explanations 
for this discrepancy: 1) The OM and UVOT UV 
filter sets have different filter transmission, 2) OM 
suffers from significantly higher level of scattered 
light than the UVOT, and third the extraction ra- 
dius of the automated OM software is 12 which is 
too large for an accurate analysis of the UV data 
due to the bright star (see UVOT section). The 
brightening of the afterglow in the last OM UVM2 
observation at 101 ks after the burst is most likely 
due to a bad subtraction of the background within 
the automatic OM data reduction software. 
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3.5. Spectral Energy Distribution 

As shown in section[331 the long-term light 
curves in all 6 UVOT filters and in X-rays do fol- 
low the same decay slope with similar break times 
at about 50 ks after the burst. In order to check 
this we determined spectral energy distributions 
(SEDs) of the afterglow at 800s, 20 ks, 100 ks, 
and 500 ks, v^rith exposure times in the XRT 
of 400s, 1.8 ks, 3.3 ks, and 4.8 ks, respec- 
tively. FigurefTOl displavs these four SEDs. These 
times are also marked in Figure[8l These times 
were picked to represent the SEDs of the earliest 
and latest time possible when the afterglow was 
detected in all 6 UVOT filters and shortly before 
and after the break. All fiuxes in all 6 UVOT 
filters and the XRT were calculated according to 
the light curves. There seem to be no obvious 
changes in the SEDs, besides the changes in the 
fluxes, over time. Another measure of any changes 
in the SEDs over time is the optical/UV to X- 
ray spectral slope or X-ray loudnes^l /3ox- For 
the afterglow of GRB 060729 we measured rest- 
frame /3ox of 0.85±0.10 at 800s, 0.84±0.05 at 20 
ks, 0.74±0.07 at 100 ks, and 0.77±0.10 at 500 ks 
after the burst. Within the errors these values 
are consistent and do not suggest any changes be- 
tween the optical/UV and X-ray parts of the SED 
over time after the first orbit. However, note that 
during the first w400s of data the SED changes 
dramatically, because the X-ray flux decays very 
fast (aa;_i=5.1 while the White and V data sug- 
gest that the optical afterglow is constant. A sin- 
gle power law spectrum between optical/UV and 
X-ray energies a /3ox=l-l is expected according to 
the fits to the X-ray data (Table[3]). This assump- 
tion of a single power law spectrum between the 
optical and X-rays is justified given that the op- 
tical has the same temporal behavior as the X- 
rays and that the optical and X-rays are both 
above the cooling frequency. The difference be- 
tween the expected /3ox=l-l and the measured 
/3ox~0.8 values suggests intrinsic reddening at the 
location of the afterglow. Based on the absorp- 
tion corrected rest-frame 2 keV flux density we 
can calculate the expected flux density at rest- 
frame 2500A. We calculated a reddening of 1.7 
mag at rest-frame 2500A which corresponds to an 



^The X-ray loudness is defined bv lTananbaum et al. 

as /3ox=-0.384 log(/2kev//„„„a )• 



ij^R-v=0.34 mag. Appl ying the relation given by 
Diplas fc Savage 1 l|l994[ n we calculated an intrin- 
sic column density A'^H.intr — 1.7x 10^^ cm^^. Con- 
sidering that this is a rough estimate, this absorp- 
tion column density agrees quite well with that 
measured from the XMM-Newton EPIC pn spec- 
trum, TVHantr = 0.85 ±0.02 x lO^i cm^^ (Table[3]). 

4. Discussion 

The afterglow of GRB 060729 has been detected 
in X-rays by the Swift-XRT longer than any other 
5'wi/i(-detected burst, up to 125 days after the 
burst. Finally by the end of December we had to 
give up on observing this burst by Swift because 
it became too faint to be detectable anymore in 
the XRT detector. Even though the afterglow was 
dropped from the Swift observing schedule after 
2006 December 27, we are still planning to obtain 
more observations with larger observatories such 
as Chandra, XMM-Newton, and Suzaku. 

4.1. Light curves 

The X-ray and UV/Optical light curves are re- 
markably similar. Not only that their decay slopes 
and break times are in good agreement (except for 
the B light curve), but they also seem to be syn- 
chronized during rebrightcning phases, which can 
be seen best in the UVW2 light curve. In par- 
ticular the re-brightening at about 15 ks after the 
burst clearly appears to be present in all 6 UVOT 
filters and in X-rays. Even though the break at 
around 50 ks after the burst seems to be achro- 
matic, we do not consider this to be a jet break. 
The post-break decay slope ax, 3=1-29 is to shal- 
low to be a jet break. If we interpret this slope as 
the post jet break slope, we would get p = 1.29, 
which is much fiatter than the electron distribu- 
tion index predicted by shock acceleration theory 
(usually p > 2). Furthermore, if p = 1.29, one 
would expect the spectral slope I3x = p/2 = 0.65 or 
l^x = (p— 1)/2 = 0.15 for a cooling frequency below 
or above the X-ray range respectively, which are 
obviously contrary to the observed spectral slope 
of /3x=l-l Most likely we have not seen the jet 
break in the afterglow of GRB 060729 because the 
afterglow has not been followed long enough as 



"The iDiplas fc Savage I l|l994l ') relation is: A^h = 4.93 X 
10^1 X Eb_v cm-2. 
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the st udies by I Willingale et all (|2007f) : ISato et al 
(|2007l ) suggest. 

One of our main results is that the hght curve of 
the afterglow does not yet show a jet break at 125 
days after the burst (81 days in the rest frame). 
This is the longest period a GRB afterglow was 
ever followed and detected in X-rays, except for 
GRB 030329 which was foll owed 258 days afte r 
the burst by XMM-Newton jTiengo et al.l 12004) 



Accor ding to the relation given in I Willingale et al 



(|2007[ ) we would have expected to see a jet break 



at about 5.5x10^ s after the burst. With an 
isotropic energy in the rest-frame 1 keV-10 MeV 
ii^iso = 1.6 X 10^ ^ ergs and the re l ations given by 



Sari et all (|l999l ) and iFrail et"all |200ll ). and the 
non-detection of a jet break up to 125 days after 
the burst, we derive that the opening angle of the 
jet has to be larger than 28°, assuming a parti- 
cle density n = O.lcm"'^ and an efficiency 77^=0.2. 
W ith ;gx=l-0 and ax =l-29 according to Table 2 
Zhang et al. ( 2006f ) we estimated an electron 



index slope p=2.3 for the ISM or wind case with 
V > max(z/m, Vc). 

The afterglow of GRB 060729 is not only re- 
markable for its long follow-up observations in 
X-rays but also for its relatively late break time 
between the flat decay phase and th e steepening 
phase (phases 2 and 3 according to IZhang et al 



20061 : iNousek et al.ll2006H Willingale et al.ll2007[ l at 
about 53 ks after the burst, which converts to 35 
ks in the rest frame. Typically the break between 
phases 2 and 3 occurs around 10 ks after th e burst 
(|Nousek et aLll2006l : Iwiihngale et al.ll2007t l. The 
late time break at 35 ks after the burst (rest frame) 
requires a substantial ongoing injection of energy 
into the afterglow. In a matter of fact, when we 
observed this afterglow with Swift e at first did not 
detect a break in the light curve until about 2 days 
after the burst due to t he lumpiness of this pla teau 
phase. Also note that Willingale et al. ( 200lf ) hst 
the break time in the X-ray light curve at 120 ks 
after the burst. 

The relatively long time of this very flat decay 
{oix,2 = 0.14 ± 0.02) phase imphes larger energy 
injection during the early time in this burst than 
in other bursts. This energy injection could be 



the result of a refreshed shock (Rees fc Meszaros 



19981 ) or continuous energy input from the cen- 
tral engine. Assume the energy injection has the 
form L{t) oc t"'', the light curve decay rate is 



(2p-4) + (p+2), ^ ^^^(^ 

20061 ). So we get q 



Vc) (jZhang et al 
0, and the total energy is in- 
creased by a factor of {Tbreak,2 /Threak,lY~'^ ~ 100 

during this energy injection phase, where Tbreak,i 
and Tbreak,2 are the first and second break time 
of the x-ray light curve respectively. Such a large 
energy increase factor is the highest amo ng the 
5w/f-detected GRBs (|Nousek et al.ll2006l ). This 
may be part of the reason that we see a bright X- 
ray afterglow for a very long time, besides that 
no jet break occurs before 125 days after the 
burst. The plateau phase of the X-ray afterglow of 
GRB 060729 is one of the longest ever observed by 
Swift. The total fluence in the 0.3-10.0 keV band 
of this plateau is about lxlO~^ ergs cm~^, which 
is about 1/3 of the total 15-150 keV fluence of the 
prompt emission. The q value {q = 0) inferred 
for this burst implies a pulsar type (i.e. magnetic 
dipole radiation) energy injectio n (e.g. iDai fc Lu 
19981 : IZhang fc Meszaros1l200l[ l 

Using the X-ray luminosity L^; = 2 x lO'^^ergs 
s~^ at t = 10 h, we estimate the isotropic ki- 
netic energy at t his time is Ek.iso ~ 2.5 x 



10^^(ej0.1)-^ergs (iFreedman fc Waxman 1 12001 



Zhang B. etall l2007f l7 where Ce is the equiparti- 
tion factor of electrons in afterglow shocks. As- 
suming the energy increasing factor 100 during 
the flat decay phase, the isotropic kinetic en- 
ergy before the energy injection phase is only 
~ 2.5 x 10^^(ee/0.1)~^ergs, implying a high ef- 
ficiency of gamma-ray production during the 
prompt phase. Using the isotropic kinetic en- 
ergy and the jet break time larger than 125 
days after the burst, we get the jet opening an- 
gle larger than Oj > 28°n]/^ (|Frail et al.ll200l[ ). 
where n_i = (n/0.1cm~'^) is the number density 
of the circumburst ISM. From this, we further 
get the beam-corrected kinetic energy of the burst 
Ekj > 1.7x 1052gj.gg(£^/o.i) 



^n]/^. Such a larger 



kinetic energy than those in usual bursts may be 
a direct consequence of the unusually long energy 
injection phase in the early time. 

Even though the X-ray and UV/optical after- 
glow of GRB 060729 was unusually bright, it was 
rather unspectacular in the BAT 15-150 keV en- 
ergy range The 15-150 keV fl uence of 2.7x10"^ 
ergs cm~^ ( Parsons et al.|[2006l ) is rather moderate 
compared to other Swift-discoYered bursts. Also 
note, that the peak luminosity of GRB 060729 of 
about 3x10^*^ ergs s~^ is very low for a long burst 
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regarding the time lag betwe en the 50-100 keV an d 
15-25 keV band as shown bv lGehrels et alJ (|2006l ). 

The observations of the X-ray and UV fields of 
GRB 060729 have been the deepest ever performed 
by Swift. In X-rays we observed the field for 1.13 
Ms and for 550 ks each in UVWl and UVW2. The 
UVWl and UVW2 observations are the longest 
exposure that were taken of any field in the UV by 
any UV observatory. The results of this study and 
the source identifications based on their spectral 
energy distributions will be presented in a separate 
paper which is in preparation. 

4.2. Spectra analysis of the early time 
data 

During the XRT WT mode epoch of obser- 
vations, two fiares are detected. For the first 
X-ray flare, only the decay part is seen by the 
XRT. The decay rates after the peak of the fiares 
are as steep as pointing to internal cen- 

tral en gine activity as the o r igin for the X-ra' 
fiares (iBurrows et al.l [2005bl: iFan fc Wei I l200. 



200i 



200i 



I 



Zhang et al.ll2006l: lOai et al.H2006l: I Falcone eTai 



Wu et a l.'2006':' Wang et al.ll2006l : lLazzati fc Perna 
Gao fc F an 200^^ 



During 130-160 s (time bins from 1 to 11 in 
Table[2|) and 190-300 s, the X-ray emission is un- 
dergoing a steep decay, which may result from the 
high-latitude emission of the corresponding fiare s 
(jKumar fc Panaitescu I l2000t iLiang et all l2006f) . 
i.e. we see the curvature effect of the radia- 
tion pulse. In this picture, the count rate and 
the peak energy (kT or -Ebrcak) both decrease 
as the pulse decays. This is because less and 
less Doppler-boosted radiation is seen from the 
pulse. This accounts for the strong correlation 
between the count rate and the blackbody tem- 
perature kT (Spearman rank order correlation co- 
efficient i?s=0.96 with Ts=11.8 and a probability 
of P< 10~^ of a random distribution; linear corre- 
lation coefficient ri=0.95). 

The WT data of the early afterglow observa- 
tions showed a dramatic change in the X-ray spec- 
trum within less than two minutes in the rest 
frame. Typically the spectra of the bins of the 
WT mode data can be fitted by a single ab- 
sorbed power law or a black body plus power law 
model. The power law model shows that there 
is a decrease in the absorption column density by 



a factor of 4 from the beginning of the observa- 
tion at 131 s after the burst to 300s after the 
burst. Such decreases in the absorption column 
density have been observed be fore in GRB after 



glows like e.g. in GRB 05071 2 (|De Pasauale et al 



I2OO6I : fLazzati fc Perna 2002), but have been usu 



ally linked to a flattening of the X-ray spectral 
slope. This is typically an artifact of the spec- 
tral fitting routine, that may be due to the cor- 
relation between spectral parameters such as A^h 
and Px- However, the situation in GRB 060729 
is completely different. Here we observe not only 
a decreasing absorption column density Nh, but 
also a steepening of the X-ray spectral slope f3x, 
which is the opposite of what one expects if this 
is just an artifact of the fitting routine. Therefore 
we consider the decrease of the absorption column 
density to be real. There are two explanations for 
a decrease in the column density of a neutral ab- 
sorber: 1) an expanding medium which results in 
a lower volume and column density and 2) ioniza- 
tion of the neutral gas. Both, the expansion and 
the ionization of the gas, happen after the explo- 
sion of the star. A softening of the X-ray spectrum 
during the initial dec line has been cominonly o b- 
served as reported by IZhang B.-B. et all (|2007t ). 

Alternatively the WT mode spectra can also be 
fitted with a black body plus power law model. In 
order to limit the number of free parameters the 
absorption column density parameter was set to 
the Galactic value. These fits show a decrease of 
the black body temperature from about 0.6 keV 
to about 0.1 keV from the beginning to the end 
of the WT observing period. Fixing the absorp- 
tion column density to a value of 1 x 10^^ cm~^, 
the value obtained at later times during the XMM 
observation, results in similar values for the black- 
body temperature. From these fits the tempera- 
ture tends to be slightly lower than when fixing 
the Nh to the Galactic value. However, within the 
errors the results are consistent. The biggest influ- 
ence the increase in the absorption column density 
has is on the normalization of the blackbody and 
power law components. 

In the latter scenario, the thermal component 
is likely to be the photospheric emission from X- 
ray fiares. In the prompt GRB phase, we have 
seen the th e rmal emission from some bursts (e.g. 



Rvde et al.l I2007I ). which has been interpreted 



as the photospheric emission when the fireball 
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becomes optically - thin (IMeszaros fc Rees 



2005 



Thompson et alJ 
Pe'er et alJl200fih . Since X- 



200C; 



2007 



Rccs & Mcszaros, 
Ramircz-Ruiz 2005 
ray flares are believed to be the result of late- 
time central engine activity with radiation mech- 
anism similar to the prompt phase, a photosphere 
component found in X-ray flares is reasonable. 
The thermal emission component discovered from 
GRB060729 also supports the internal origin of 
the X-ray flares, rather than external shocks. Let 
us examine the relation between the black body 
radius Rbb derived from the spectral fitting and the 
photosphere radius Rph ■ For a relativistic moving 
source, the luminosity of the thermal component 
at the photospheric radius is 



L = 4TTRl^r^aT"^ 



(1) 



where T' is the photospheric temperature in the 
co-moving frame and T is the bulk Lorcntz fac- 
tor. As the usual fitting uses L = AirRlf^aT'^, we 
get Rph = TRbb by taking advantage of T = FT'. 
The fitted black body radii around the peak of the 
X-ray flares are a few times lO^^cm. This means 
the photospheric radii of the two X-ray flares in 
GRB060729 are a few times of lO^^Fcm. The 
photospheric radius is very se nsitive to the bulk 



Lorentz factor F of the fireball ( Rccs fc Mcszaros 



20051 ) and a photospheric radius of 10 
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lOi^c 



is reasonable if F of the X-ray flare is of the order 
of ten. 

During the early steep decay phase, the bolo- 
metric flux (equivalent to the flux in XRT band 
if the peak energy located within the 0.3-10 kcV 
range of the XRT) may decrease as Fx{t) oc 
(t — to)~^ due to the curvature effect (Rydc & 
Petrosian 2002), where to is some reference time of 
the flare (Liang ct al. 2006), and the peak energy 
decays as (t — io)^^- So if we fit the spectrum 
with the black body model {Fx{t) oc T^Rl^) all 
the time during the decay phase, we would expect 
that the black body radius increases with time as 
Rhb oc {t — to). This may explain the apparent in- 
crease of the black body radius with time during 
the steep phase of the flares. 

5. Conclusions and Summary 

We studied the Swift and XMM-Newton obser- 
vations of the afterglow of GRB 060729 and found: 

• The light curve of the afterglow extends out 



to 125 days after the burst (81 days in the 
rest-frame) without showing any sign of a jet 
break. We estimated that the jet opening 
angle has to be larger than 28°. This is the 
longest foUowup and detection of a GRB af- 
terglow in X-rays ever performed, except for 
GRB 030329. 

• The X-ray light curve can be generally 
described by an initial steep decay slope 
ai=5.1±0.2, a break time Tbroak.i=530±25 
s, flat decay slope a2=0.14±0.02 with a 
break at Tbrcak,2=56.8±10 ks, and a steep 
decay slope a3 = 1.29±0.03. 

• The unusually long flat decay phase of the af- 
terglow of GRB 060729 implies a much larger 
energy injection than seen in any other GRB 
afterglow. 

• After the initial phase, the light curves in X- 
rays as well as in all 6 UVOT filters follow 
the same shape. 

• In the initial phase the afterglow shows a 
dramatic change in its X-ray spectrum which 
can either be described by a steepening of 
a power law spectrum with a simultaneous 
decrease in the intrinsic column density, or a 
decrease in the black body temperature from 
0.6 keV at 130s after the burst to 0.1 kcV at 
250s observed after the burst. 

• The spectral analysis of the Swift XRT PC 
mode and XMM-Newton EPIC pn and MOS 
data shows that the X-ray spectrum of the 
afterglow agrees with an absorbed power law 
with /3x=l-l and an intrinsic column density 

iVH,intr = 1 X 1021 Cm'^. 

• The reddening and intrinsic column density 
estimated from the spectral energy distribu- 
tion agrees well with the value found from 
the XMM-Newton analysis. 
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Fig. 1.— UVWl image (5f7 x 4.'2) of the field of GRB 060729 with an exposure of 550 ks. The circles at the 
source position displays the XRT position as given in Section lXTl The circles in the large image show a 8 
radius at the source just for display purposes, and the 8 background extraction region. The zoom-in image 
(80 X 60 ) in the upper right corner shows the 3'.'5 XRT error radius. Note the bright star which is only 
107 away from the position of the afterglow of GRB 060729. 



GRB 060729 BAT Mask-weighted Light Curve 




-100 100 200 

Time from BAT trigger (sec) 



Fig. 2.— Swift BAT light curves in the 15-25 keV band (top), 25-50 keV (second panel), 50-100 keV (second 
lower panel), and in the 15-100 keV band (bottom). The vertical dashed line in the 15-100 keV plot marks 
the time for which Tgo was calculated and the dotted-dashed line the beginning of the XRT observations at 
124s after the burst. 
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Table 1 

Spectral fits to the BAT data. The BAT data have been divided into 5 segments as listed 

BELOW. 



„ 2 / Flux^ Flucncc"^ 

Lggor /3i5-i50kev X Iv i5_i5o koV 1 kcV - 1 McV^ 1 koV - 10 McV^ 15-150 kcV 1 koV - 1 MgV=' 1 koV - 10 MeV^ 



1st peak 


- 


10 


1 05+"*= 


57/56 


3.2x10"** 


1.0x10"^ 


1.3x10"'' 


3.2x10"' 


1.0x10"'' 


1.3x10"" 


2nd peak 


70 - 


- 88 


0.59l°:ii 


59/56 


5.8x10"** 


1.6x10"^ 


4.5x10"'' 


1.0x10"'' 


2.8x10"° 


8.1x10"" 


3rd peak 


88 - 


124 




50/56 


2.8x10"** 


7.7x10"** 


1.2x10"'' 


1.0x10"" 


2.8x10"" 


4.3x10"" 


XRT flare 1 


124 - 


160 


-^■'^"-0.86 


49/56 


2.6x10"" 


1.1x10"** 


1.3x10"** 


9.4xlO~* 


4.0x10"'' 


4.7x10"'' 


XRT flare 2 


180 - 


190 




53/56 


1.8x10"'^ 


1.5x10"** 


1.5x10"** 


1.8x10"** 


1.5x10"' 


1.5x10"' 



15-150 kcV fiux in units of ergs s cm ^ 
"^Fluencc in units of ergs cm""^ 

^Rest-frame 1 kcV - 1 MeV (0.65 keV - 650 keV observed) and 1 kcV - 10 MeV (0.65 keV - 6.5 MeV observed) 



Table 2 

Spectral analysis of the 21 bins of the Swift WT mode data.^ 















Sit 


iglc Power 


Law 






Blackbodv + Power Law 




Power Law with Exponential Cutoff 


Bin # 


Time 




















bb-Flux^ 




X = /^ 


^cut of f ^ 




2 


132 


1.2 


919ib30 
841±28 


55zb0 06 
0.60±0.06 


4, 


'^•^-1.15 


1.51+;; 


.31 


21/23 





=■"-0.04 


7.46 


2.49 


16/23 


n so+0.44 
O'^^-O.OS 


15/23 


3 
4 


133 
135 


1.4 
1.6 


807±27 
737±25 


0.51±0.05 
0.43zb0.05 


4. 


.1+0.88 
■^-^-0.77 


2.10+« 


29 
.25 


39/32 





45 + 0.04 
-0.03 


6.13 


3.46 


51/32 


-„ + 0.36 
0"=-0.02 


48/32 


5 
6 


136 
138 


1.6 
1.8 


712±24 
635±22 


0.47±0.05 
0.42±0.05 


5, 


— 0.85 


2.40+0 

— 


33 
.30 


36/31 





,43 + 0-02 
— 0.02 


6.23 


3.91 


34/31 


0.51+9-=?! 
— 0.39 


33/31 


7 

8 


140 
142 


2.1 
2.6 


579±19 
520±18 


0.33±0.04 
0.29±0.04 


4. 


..^ + 0.58 
^■^-0.53 


2.47^0 


24 
.21 


39/46 





+ 0.02 
■^^-0.02 


5.72 


4.69 


45/46 


+n 46 

u.oj_o 25 


45/46 


9 
10 


145 
149 


3.1 
4.4 


498±16 
380±12 


0.25±0.04 
0.12±0.04 


4. 


70+0.65 
"'-0.58 


2.99+° 


32 
.28 


75/50 





+ 0.02 
■^^-0.02 


3.60 


5.48 


96/50 


+ 3 
0.50102 


107/50 


11 


154 


5.7 


289±10 


-0.00±0.04 


3. 


""-0.56 


2.64+0 


34 
.30 


29/34 





+ 0.02 
'^^-0.02 


2.81 


7.35 


30/34 


67+0-*'' 

"■"'-o.og 


33/34 


12 


161 


7.3 


227±7 


-0.07±Q.04 


3. 


" -0.63 


2.73 + 


41 
.36 


43/32 





29+0.02 
'^''-0.02 


1.98 


4.84 


32/32 


64+0-'"' 


34/32 


13 


167 


6.4 


258±9 


0.10±0.04 


2, 


,„+0.60 
^■'-0.54 


2.04+0 


33 
.29 


38/32 





,-.+ 0.03 
-0.03 


1.45 


4.77 


43/32 


1 15+0-61 
^--'^-0.34 


48/32 


14 


173 


5.2 


320±10 


0.21±0.04 


3. 


+ 0.62 
°"-0.56 


2.46+0 


27 
.24 


40/33 





34+0.03 
^ -0.02 


3.54 


4.68 


39/33 


72+0-37 


37/33 


15 


179 


5.4 


310±10 


0.15±0.04 


2, 


., + 0.59 
*^-0.53 


2.01 + 


30 
.28 


39/35 





o,+0-03 
-0.02 


2.90 


4.61 


37/35 


-O+0.40 
0-™-0.03 


35/35 


16 


184 


5.5 


301±10 


0.02±0.04 


2, 


97+0.67 
^ -0.60 


2.571S 


38 
.34 


35/34 





,„+0.02 
"^"-0.02 


3.28 


7.52 


32/34 


^„+0.56 
0-58-0.23 


36/34 


17 


190 


7.4 


227±7 


-0.09±0.03 


3. 


+ 0.85 
*''-0.76 


3.07j:0 


51 
.46 


55/35 





+ 0.01 
-^^-0.01 


3.24 


8.76 


40/35 


51 + 0-18 


45/35 


18 


199 


9.4 


178±5 


-0.26±0.03 


2, 


4g + 0.62 
*°-0.55 


2.98 + 


41 
.36 


62/39 





19+0-01 
^"-0.01 


3.44 


14.84 


58/39 


,5+0.62 
0-5°-0.29 


60/39 


19 


210 


14.0 


127±4 


-0.31±Q.03 


2, 


,,+0.63 


3.07±0 


42 
.36 


59/38 





+ 0-01 
^»-0.01 


2.84 


15.40 


33/38 


51+0-69 

"-'"-o.si 


60/38 


20 


229 


23.4 


72±2 


-0.50±0.03 


1, 


„, + 0.53 
= ^-0.44 


3.28 + 


41 
.35 


72/36 





+ 0.01 
^-0.01 


2.15 


18.19 


38/36 


51 + 0-6 


76/36 


21 


296 


112.0 


14.8±0.4 


-0.56±0.03 


0. 


g^+0.35 
° -0.30 


2.83 + 


36 
.32 


64/37 





,,+0.01 
^^-0.01 


0.50 


15.65 


52/37 







"'"For the blackbody plus power law and the power law with exponential cutoff model the power law slope 
parameter was fixed to the Galactic value [4.82 X 10^0 cm"^ [Jjickey fc Loekm'arnil99(ji ). 

^Count rate CR in units of counts s'"'"'". 

^The hardness ratio is defined as H R = (hard — so ft) / (hard -|- soft) where soft and hard are the ph 
^The column density iVjj is given in units of 10'^"'" em~~^. 
^Blackbody temperature in units of keV. 

^The fluxes are given in units of 10""^ ergs s""'" cm~'^. 



! fixed to /3^ = 1.0. For the blackbody plus power law model the absorptic 



Blackbody radius J?]^}-, given in units of 10 
^The break energy is given in units of keV. 



12 , 



, the 0.3-1.0 and 1.0-10.0 keV band. 




Time since burst [s] 

Fig. 3.— Combined Swift BAT and XRT WT light curves. The Figure clearly shows that the XRT was 
starting observing GRB 060729 during the fourth peak seen in the BAT. 



Table 3 

Spectral fits to the Swift XRT PC mode and XMM EPIC pn and MOS data 



Detector 


7Vh = 


/3x 








2 / 


Swift XRT (20-40 ks after burst) 

Swift XRT (44.9-107.5 ks after burst) 

Swift XRT (200 ks after burst) 

XMM EPIC pn 

XMM MOS1+MOS2 

XMM MOS 1+2 + Swift XRT 

XMM pn + MOSl+2 + Swift XRT^ 


14.65t«:- 
8.58;S:f„ 

9 54+°-^= 

°->"-0.17 


, „,+0.10 

1 lQ+0.09 

, ,7+0.25 
^■^ ' -0.22 
, ,,+0.01 
^■^^-0.01 

-^■"^-0.02 
, r,t;+0.02 
1-U3_Q 02 

1 ns+"''i 
1.08_Q 01 


87/103 
134/120 
33/25 

1159/1149 
937/800 
1106/920 

2599/2071 


21 57+*-22 
19.20+^-™ 
17 46+ " 

7 70 + 0.39 
' ■ ' ^-0.38 
Q oQ+0.59 
'^■■^'^-0.58 
S K1+0.06 

7 cr,+0.03 
' -^^-0.03 


1.12±0.08 

1 1 1 +0.08 
-^■-^-^-0.07 

-^■-^"-0.19 

-^■-^-'--o.oi 

^•"^-0.01 
^•'^''-0.01 


88/103 
134/120 
33/25 
1124/1149 
929/800 
1097/920 
2508/2071 



^XMM observed the afterglow between 44.9-107.5 ks after the burst. 
^The column density A^h at z— is given in units of lO'^^ cm""^. 

^Intrinsic column density A^H,intr at the rcdshift of the burst, z— 0.54, is given in units of lO^'^ cm^'^. The absorption 
column density at z— is set to the Galactic value, 4.82x10'^'^ cm~'^. 
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Fig. 4. — Swift XRT light curve of the Windowed Timing (triangles) and photon counting (crosses). The 
downward arrow marks the Scr upper limit at the end of the Swift observations. This upper limit contains 
a total exposure time of 63.5 ks obtained between 2006 December 8 and December 27. The dotted vertical 
lines mark the start and end times of the XMM-Newton observation. 
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Fig. 5. — Swift XRT Windowed Timing mode light curve. The panels display from the top to the bottom 
the XRT count rate (in units of counts s""^), the hardness ratio (see text for definition), the X-ray spectral 
slope /3x of a single power law fit, the free fit column density Nh in units of 10^^ cm~^, the blackbody 
temperature kT, the blackbody radius i?bb (hi units of 10^^^ cm), and the cut off energy E'brcak of a power 
law with exponential cutoff. All these fit parameters are listed in Table [2l The numbers in the top panel 
mark the bins that were used for the spectra shown in Figure [6l 
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0.5 1 2 5 

Energy [keV] 



Fig. 6. — Swift XRT WT spectra of bins 1+2 (black filled circles), 12 (red open squares), 14 (green open 
circles), and 21 (blue filled squares) fitted by absorbed single power laws as given in Table[2] The numbers 
refer to the bins as shown in Figure El 
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■". _. . . . . . I . 



Fig. 7. — Absorbed power law model fits with 7Vh=8.57x 10^° cm"^ and /3x=1.08, as listed in Tabled to 
the XMM-Newton EPIC pn (black filled circles), MOSl (red open squares), M0S2 (green open circles), and 
Swift XRT PC niode(blue filled squares) spectra. The Swift XRT data were selected between 44.9 ks to 107.5 
ks after the burst, simultaneous with the XMM-Newton observation. 
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Fig. 8.—- Swift UVOT and XRT PC mode light curves, with the XRT light curve on the top and the UVOT 
light curves starting with UVW2 on the top and V on the bottom. The downward arrows mark the 3a upper 
limits as listed in Tabled The vertical lines mark the times when the spectral energy distributions of the 
afterglow were determined as shown in Figure fTOl 
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Table 4 

Central times, exposure times and aperture corrected magnitudes of the UVOT light 

CURVES 



UVWl-Filtc 



UVM2-Filtc 



UVW2-FiltG 



# 


Timol 


2 

cxp 




Mag 




2 

cxp 


ag 


Timol 


-'cxp 




Mag 


Timcl 


2 

cxp 




M 




Timcl 






Mag 


Timcl 


2 

cxp 


1 


120 


9 


17 


07±0.31 


715 


10 


18.27±0.38 


696 


20 


16 


48±0.12 


673 


20 


16 


92±0 


20 


649 


18 


16 


Sl±0.24 


749 


20 


2 


440 


390 


17 


37±0.06 


1451 


20 


18.03±0.23 


844 


20 


16 


88±0.15 


820 


20 


16 


S9±0 


20 


796 


20 


17 


57±0.33 


1489 


20 


3 


773 


19 


17 


36zb0.26 


1609 


20 


18.10±0.26 


1427 


20 


16 


85±0.16 


1403 


20 


17 


47±0 


26 


1379 


20 


17 


12±0.26 


1647 


20 


4 


1164 


393 


17 


32±0.06 


1767 


20 


18.26±0.30 


1585 


20 


16 


85±0.15 


1561 


20 


17 


31dzO 


24 


1537 


20 


16 


85±0.23 


1805 


19 


5 


1513 


19 


17 


45±0.32 


6224 


197 


18.45±0.10 


1743 


20 


16 


49±0.13 


1719 


20 


17 


14±0 


23 


1695 


20 


16 


75±0.22 


6634 


197 


6 


1671 


19 


17 


16±0.24 


7657 


197 


18.45±0.11 


6020 


197 


17 


52±0.07 


5815 


197 


17 


55±0 


09 


1850 


13 


17 


63±0.45 


7998 


63 




1829 


19 


17 


55±0.34 


12919 


134 


18.68±0.14 


7452 


197 


17 


54±0.07 


7247 


197 


17 


62±0 


09 


7043 


197 


17 


80±0.12 


11878 


751 


8 


6838 


197 


18 


10±0.16 


18042 


211 


17.86±0.06 


12778 


134 


17 


89±0.11 


12568 


268 


18 


09±0 


10 


13880 


377 


17 


76±0.09 


13266 


537 


9 


13613 


134 


18 


38±0.25 


23824 


211 


17.84±0.06 


17822 


211 


16 


91±0.05 


17495 


422 


17 


11±0 


05 


19550 


600 


17 


O7±0.O5 


18585 


844 


10 


19128 


211 


17 


33±0.09 


29606 


211 


17.78±0.06 


23603 


211 


17 


20±0.06 


23276 


422 


17 


23±0 


05 


25338 


599 


17 


20±0.05 


24367 


845 


11 


24911 


211 


17 


73±0.11 


35388 


211 


18.14±0.08 


29387 


211 


17 


17±0.06 


29059 


422 


17 


22±0 


05 


31109 


600 


17 


37±0.O6 


30149 


844 


12 


30692 


211 


17 


53±0.10 


41205 


208 


18.12±0.08 


35167 


211 


17 


20±0.06 


34840 


422 


17 


35±0 


05 


36892 


599 


17 


40±0.06 


35931 


845 


13 


36475 


211 


17 


98±0.13 


47494 


147 


18.22±0.09 


40989 


207 


17 


26±0.06 


40666 


415 


17 


35±0 


05 


42681 


587 


17 


45±0.06 


41738 


829 


14 


42273 


207 


17 


86±0.13 


56873 


129 


18.31±0.11 


47340 


147 


17 


18±0.07 


47109 


296 


17 


24±0 


06 


48550 


416 


17 


25±0.06 


47877 


592 


15 


48259 


147 


17 


47±0.11 


67943 


42 


18.94±0.33 


53746 


83 


17 


32±0.10 


52836 


183 


17 


36±0 


08 


54427 


229 


17 


63±0.10 


54050 


330 


16 


54265 


83 


18 


12±0.26 


75864 


211 


18.77±0.12 


62872 


69 


17 


24±0.11 


58802 


385 


17 


60±0 


06 


63155 


177 


17 


93±0.14 


62989 


277 


17 


63099 


69 


18 


29±0.35 


81646 


211 


18.88±0.13 


69860 


211 


17 


75±0.0S 


64698 


436 


17 


76±0 


07 


100500 


603 


18 


66±0.11 


76105 


249 


18 


91349 


223 


18 


68±0.22 


87426 


211 


18.96±0.14 


75644 


211 


17 


84±0.0S 


69533 


422 


17 


85±0 


07 


106280 


599 


18 


57±0.11 


82086 


640 


19 


100080 


211 


18 


45±0.22 


93208 


211 


18.72±0.11 


81426 


211 


18 


07±0.10 


75315 


423 


17 


93±0 


07 


112060 


597 


18 


58±0.11 


87970 


846 


20 


105860 


211 


18 


87±0.28 


98988 


211 


19.11±0.15 


87206 


211 


17 


98±0.10 


81097 


423 


18 


05±0 


08 


117840 


600 


18 


52±0.10 


93752 


846 


21 


111650 


211 


18 


49±0.19 


104770 


211 


19.02±0.13 


92987 


211 


18 


03±0.10 


86877 


423 


18 


17±0 


08 


123630 


600 


18 


90±0.13 


99532 


847 


22 


117430 


211 


19 


03±0.32 


110560 


211 


18.84±0.13 


98767 


211 


18 


40±0.12 


92658 


423 


18 


35±0 


09 


129410 


601 


18 


81±0.12 


105320 


846 


23 


123210 


211 


18 


40±0.20 


116340 


211 


19.18±0.16 


104550 


211 


18 


23±0.11 


98438 


423 


18 


33±0 


09 


135340 


294 


18 


79±0.17 


111100 


846 


24 


128990 


211 


18 


92±0.27 


122120 


211 


19.62±0.23 


110340 


211 


18 


50±0.14 


104230 


423 


18 


53±0 


10 


146960 


342 


18 


99±0.19 


116880 


846 


25 


135120 


111 


18 


77±0.38 


127900 


211 


19.32±0.17 


116120 


211 


18 


33±0.12 


110010 


423 


18 


52±0 


10 


195280 


737 


19 


62±0.20 


122660 


845 


26 


158090 


156 


19 


36±0.56 


140600 


225 


19.37±0.19 


121900 


211 


18 


60±0.16 


115790 


423 


18 


56±0 


10 


204280 


189 


19 


48±0.11 


128450 


845 


27 


226950 


941 


19 


82±0.29 


154780 


180 


19.77±0.30 


127680 


211 


18 


43±0.12 


121570 


422 


18 


63±0 


11 


235840 


471 


19 


84±0.29 


134830 


445 


28 


319820 


1374 


20 


22±0.34 


168380 


211 


19.66±0.24 


134430 


110 


18 


46±0.17 


127360 


423 


18 


76±0 


11 


288290 


1550 


19 


86±0.16 


146750 


522 


29 


576790 


5531 


21 


04±0.40 


174170 


209 


20.00±0.32 


146620 


178 


19 


08±0.25 


137350 


412 


18 


79±0 


12 


331760 


2240 


20 


21±0.18 


168860 


729 


30 


897975 


2368 


3cr 


ul=20.92 


179940 


211 


20.00±0.32 


156590 


211 


18 


64±0.15 


145550 


291 


18 


76±0 


14 


397970 


3140 


20 


33±0.16 


174710 


839 


31 










194290 


275 


19.68±0.20 


168160 


211 


18 


70±0.15 


150480 


423 


18 


78±0 


12 


487880 


2980 


20 


59±0.20 


180340 


556 
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206460 


273 


19.29±0.15 


173950 


209 


18 


84±0.16 


154560 


456 


18 


94±0 


12 


663250 


8830 


21 


68±0.30 


194650 


1100 


33 










319120 


1374 


20.11±0.14 


179720 


211 


19 


09±0.21 


162860 


261 


19 


18±0 


19 


926689 


3976 


3a- 


ul=21.96 


206810 


1092 


34 










449360 


2304 


20.38±0.13 


194150 


275 


18 


82±0.13 


167830 


423 


19 


04±0 


13 
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662500 


3619 
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206310 
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19 


03±0.16 


173630 


420 
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98±0 


13 










250410 


1111 
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897470 


2368 
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221160 
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2304 
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220980 
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15 
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286960 


1100 
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Fig. 9.— Swift UVOT White and V and XRT WT PC mode early time light curves of the first orbit. The 
left panel shows the UVOT White and the XRT WT mode data, and the right panel shows UVOT V and 
XRT PC mode data. The binning in the PC mode data is 25 photons per bin. The UVOT White data are 
binned in 10 s intervals and the V data are binned in 25 s intervals. 



Table 5 

Decay slopes^ and break times of the UVOT light curves. 



Filter 02^ Tb^oak^ "3 



V 


0.24±0.05 


40+15 


1.21±0.09 


B 


0.16±0.05 


,(,+20 


0.96±0.05'' 


U 


0.41±0.07 


54+10 


1.40±0.07 


UVWl 


0.47±0.06 


48+5 


1.29±0.03 


UVM2 


0.35±0.10 




1.39±0.05 


UVW2 


0.46±0.07 


00 -10 


1.36±0.04 


X-rays 


0.14±0.02 


57 + 10 


1.29±0.03 



^ The decay slop es 02 and follow 
the conventions in iZhang et al.| i2006r and 
INousek etim | |2006D . 

^Decay slope before the energy injection and 
rebrightening at 20 ks after the burst 

^Time after the break in ks in the observed 
frame 

"^Note that when limiting the analysis of the 
late time decay slope to 50-200 ks after the 
burst, the decay slope is 1.17±0.16. 
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Fig. 10. — Spectral energy distributions of the afterglow of GRB 060729 at 800s (upper left), 20 ks (lower 
left), looks (upper right) , and 500 ks (lower right). The UVOT photometry data are corrected for Galactic 
reddening (.ER_v=0.050 lSchlegel et al.lll998l ). 
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Table 6 

List of the Swift UVOT White and V observations during the first orbit as shown in 

FigureIH 



Time^ 


White 


V 


140 


1 '7 T7_Ln 1 O 

1 1 .2 1 ±U. io 




150 


1? F;i_|_n 11 
i 1 .QlltU.Zi 




160 


1 '7 n o _L n oT 
i ( .yozhU-Z ( 




170 






180 


io.UiitU.oU 




190 


ocrul: lo.o/ 




200 


ocrui, io.ou 




210 






220 


io. i / itU.oO 




230 


i f .44±U. i 




278 




1 / .OUihU.oO 


303 




Q , 1 , If? O 1 


328 




i / .OZItU.ZD 


353 




q , 1 , 1 frn 
OCTUl, 1 I . iu 


378 




Q , 1 , 1 ? fTQ 

OCTUl, 1 1 . li) 


403 




17 1 'T-l-n 0/1 
i i . 1 1 ItU.Z4 


428 




ic no_l_n on 


453 




17.03±0.22 


478 




17.26±0.27 


503 




17.12±0.23 


528 




16.88±0.19 


553 




17.42±0.31 


578 




17.20±0.26 


603 




17.46±0.30 


628 




16.89±0.21 


778 




17.27±0.28 


978 




17.59±0.34 


1003 




17.66±0.33 


1028 




17.55±0.35 


1053 




16.94±0.19 


1078 




17.60±0.32 


1103 




17.02±0.21 


1128 




17.40±0.29 


1153 




3ctu1: 17.73 


1178 




17.37±0.27 


1203 




16.97±0.20 


1228 




17.43±0.29 


1253 




16.86±0.19 


1278 




3aul: 17.70 


1303 




16.94±0.19 


1328 




17.42±0.29 


1353 




17.08±0.22 


1503 




3(Tul: 17.19 


1528 




3ctu1: 17.81 


1678 




16.96±0.27 



^The times note the middle of the ob- 
servation in s after the burst 
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Table 7 

List of the XMM-Newton Optical Monitor observations 



ObsID Tobs^ Filter Mag 



006 


48351 


U 


17.27±0.02 


018 


52707 


U 


17.43±0.02 


010 


70984 


UVWl 


17.53±0.04 


015 


75291 


UVWl 


17.57±0.04 


016 


84105 


UVWl 


17.71±0.04 


020 


88612 


UVWl 


17.83±0.05 


012 


92919 


UVM2 


17.89±0.12 


014 


97226 


UVM2 


17.98±0.13 


017 


101533 


UVM2 


17.73±0.11 



^The times note the middle of the observa- 
tion in s after the burst 
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